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1.  INTRODUCTION 


1 . 1 Background 

Vessels  containing  gases  at  high  pressures  are  used  throughout  industrial  and 
research  areas.  Their  users  range  from  gasoline  service  stations  to  nuclear  power 
plants  to  spacecraft.  The  potential  hazards  from  accidental  rupture  must  be  con- 
sidered in  their  design  and  siting.  Hazards  include  blast  waves  and  high  energy 
fragments  generated  by  vessel  rupture. 


Blast  and  fragment  data  from  the  pneumatic  burst  of  high  pressure  vessels  are 
not  extensive.  Empirical  data  include  the  low  pressure  work  of  reference  (l)  in 
which  an  Atlas  missile  was  burst  at  a pressure  of  1*9-7  psiA  and  the  work  of 
reference  (2)  which  is  based  on  the  rupture  of  small  glass  spheres  at  pressures  up 
to  1*00  psiA.  The  higher  pressure  work  in  reference  (3)  is  based  on  the  burst  tests 
of  full  scale  vessels  at  pressures  ranging  from  61+0  to  8ll*5  psiA.  Data  from 
vessels  burst  above  8ll*5  psiA  are  not  available. 

The  investigation  described  in  this  report  was  designed  to  extend  empirical 
data  to  higher  pressures  and  to  generalize  the  results  through  computer  techniques. 


1.2  Scope  of  the  Investigation 


Seven  vessels  were  pressurized  with  argon  until  they  burst.  The  area  was 
instrumented  to  measure  the  blast  and  fragment  parameters  generated  by  vessel 
rupture.  All  vessels  were  T-l  steel  spheres  with  an  internal  volume  of  1 ft  . 
Design  burst  pressures  were  15,000,  30,000,  and  50,000  psi. 


Blast  calculations  were  made  in  two  ways.  A one-dimensional  calculation  was 
used  to  calculate  airblasr  parameters  in  a fiee-air  situation.  A two-dimensional 
hydrocode,  TUTTI,  was  developed  to  give  a better  calculation  of  the  blast  field 
about  vessels  burst  in  real  situations. 


One  vessel  was  burst  inside  a five-sided  steel  structure  to  obtain  blast 
loading  on  containment  walls  in  a real  situation.  The  two  dimensional  hydrocode 
can  be  used  to  obtain  blast  loads  on  a confining  structure's  walls. 


(1)  Moskowitz,  H. "Blast  Effects  Resulting  from  Fragmentation  of  an  Atlas  Missile", 
AIAA  Paper  No.  65-195,  AIAA/NASA  Flight  Testing  Conference,  Huntsville,  Alabama, 
February  15-17,  1965 

(2)  Boyer,  D.  W.  et  al  "Blast  from  a Pressurized  Sphere",  University  of  Toronto 
Institute  of  Aerophysics,  Report  No.  1*8 , January  1958 

(3)  Pittman,  J.  "Blast  and  Fragment  Hazards  From  Bursting  High  Pressure  Tanks", 
Naval  Ordnance  Laboratory,  Report  No.  N0LTR  72-102,  May  17,  1972 

Q 

UNCLASSIFIED 


- 1 


Y 


UNCLASSIFIED 

NSWC/WOL/TR  75-87 


EXPERIMENTAL  OUTLINE 


The  experimental  plan  is  given  in  detail  in  Annex  A.  It  is  summarized  here  to 
aid  in  understanding  the  investigation. 

For  the  investigation,  a pressure  vessel  was  located  a few  inches  above  ground/ 
firing  pad  level  and  was  pressurized  with  argon  until  it  burst.  The  \essel  was 
instrumented  to  monitor  internal  ga3  pressure  and  temperature  and  circumferential 
expansion.  The  area  around  the  vessel  was  instrumented  to  measure  the  airblast 
generated  by  the  rupturing  vessel.  A system  of  breakwires  was  assembled  about  the 
vessel  to  measure  vessel-wall  fragment  velocities.  A semi-circular  arena  made 
from  pressed  board  panels  was  set  up  to  recover  vessel  wall  fragments. 

The  expermental  layout  is  shown  in  Figures  2.1  and  2.2.  A drawing  of  a 
50,000-psi  vessel  is  shown  in  Figure  2.3.  Details  of  all  pressure  vessels  are 
given  in  Table  2.1. 

3.  AIRBLAST  RESULTS 

Airblast  generation  by  high  pressure  vessel  rupture  is  controlled  by  a number 
of  factors.  These  include  the  geometry  of  vessel  breakup,  the  volume  of  the 
vessel,  the  type  of  gas  in  the  vessel,  and  the  temperature  and  pressure  in  the  ^ 
vessel  at  failure.  In  this  experiment,  the  test  vessel  internal  volume  was  1 ft 
for  all  tests.  The  internal  gas  temperature  at  rupture  was  between  300  and  371*°K, 
a spread  not  expected  to  affect  airblast  generation.  A monotomic  fill  gas,  argon, 
was  used.  Thus  the  only  variables  in  the  experiment  were:  the  controlled  one, 

vessel  pressure  at  rupture,  and  the  uncontrolled  one,  the  geometry  of  vessel 
breakup. 

3.1  Effects  of  Rupture  Geometry 

All  the  vessels  used  in  this  test  program  burst  into  two  pieces.  Both  the 
15,000  and  30,000-psi  vessels  burst  into  two  hemispheres.  The  50,000-psi  vessels 
separated  into  two  pieces,  one  slightly  larger  than  the  other  as  discussed  in 
Annex  B where  vessel  rupture  histories  have  been  reconstructed.  From  these 
reconstructions,  we  determine  the  direction  of  maximum  gas  motion,  or  jetting. 
Asymmetries  in  the  blast  field  about  the  bursting  vessel  correlated  with  the 
direction  of  the  Jetting. 

The  estimated  direction  of  jetting  for  each  shot  is  plotted  in  Figure  3.1. 

This  estimates  the  Jet  direction  in  the  horizontal  plane  only.  Estimates  of 
jetting  in  the  vertical  direction  will  be  discussed  for  specific  cases. 

Airblast  overpressure  data  from  all  unconfined  tests  are  plotted  in  Figures 
C.5,  C.9,  and  C.13  in  Annex  C.  Curves  fitted  by  eye  to  these  data  are  shown  in 
Figures  3.2,  3.3,  and  3.1*.  These  figures  show  the  overpressure/ground  range 
relationship  for  gage  lines  1 and  2 for  each  test.  A single  gage  was  deployed  on 
a third  line  after  shot  1.  Data  recorded  at  this  position  were  suspect  for  shots 
3 and  L in  that  it  generally  read  half  the  expected  overpressure.  The  entire  gage/ 
electronics  were  changed  after  shot  after  this  overpressures  were  recorded  that 
support  the  conclusions  drawn  here. 
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Curves  from  Figures  3.2,  3.3,  and  3.b  have  been  replotted  in  Figures  3.3,  3.6, 
and  3.7;  this  time  they  are  grouped  according  to  their  angular  relationship  to  the 
estimated  Jet  direction.  Figure  3.5,  shows  overpressure/ground  range  curves 
measured  along  lines  ranging  from  13c  to  33°  from  the  estimated  direction  of 
maximum  argon  motion.  At  ground  ranges  greater  than  20  feet,  (overpressures  less 
than  5 psi),  the  vessel  burst  pressure  has  little  effect  on  airblast  overpressures. 

Even  at  closer  ranges  and  higher  overpressures,  differences  in  overpressure  in 
Figure  3.5  may  have  been  due  to  geometry  rather  than  burst  pressure.  For  instance, 
the  15,000-psi  vessel's  bottom  half  served  as  a shield  for  the  closer  gage  positions 
(see  Figure  B-5).  For  the  22°,  50,000-psi  curve, the  vessel  burst  so  as  to  direct 
the  blast  upward.  The  33°,  50,000-psi  vessel's  blast  was  directed  downward  (see 
Figure  B-12  and  the  paragraph  below). 

The  angles  between  the  gage  lines  and  the  argon  jets  were  greater  for  the  data 
displayed  in  Figure  3.6.  The  anomaliy  where  the  68°,  50,000-psi  curve  shows  higher 
pressures  than  the  57°,  50,000-psi  curve  may  be  due  to  the  fact  that  the  first 
opening  for  one  tank  occurred  along  a line  described  by  the  tanks  center  and  a 
point  on  its  surface  that  was  nearly  vertical  to  the  ground  in  one  case  (57°)  and 
horizontal  in  the  other  (68°).  Thus  for  the  68°  case,  close  in  gages  were  shielded 
giving  low  pressure  close  in;  but  the  Jet  was  broader  horizontally,  tending  to 
generate  higher  pressures  further  out.  The  57°  jet  was  broad  in  the  vertical 
direction.  The  family  of  curves  in  Figure  3.6  indicate  lower  overpressures  for  a 
given  distance  than  is  shown  in  Figure  3.5.  They  are  however,  higher  than  those 
for  even  larger  angles  shown  in  Figure  3.7. 

The  overpressure/ground  range  relations  plotted  in  Figure  3.7  show  a greater 
angular  distance  from  the  jet  than  is  shown  in  Figures  3.5  and  3.6;  and  therefore, 
lower  pressures  were  measured  along  these  lines  than  those  in  Figures  3.5  and  3.6. 

The  greater  angular  dispersion  of  the  data  in  this  figure  also  make  for  wider 
dispersion  of  the  curves. 

From  the  above,  the  factor  controlling  the  magnitude  of  the  airblast  overpressure 
from  the  high  pressure  vessels  ruptured  in  this  investigation  was  the  -Jetting 
direction  of  the  high  pressure  gases.  For  fill  gas  temperatures  in  the  range  of 
300°K,  increasing  rupture  pressure  does  not  significantly  increase  the  airblast 
pressure  for  rupture  pressures  above  15,000  psi. 

3.2  Effects  of  Temperature  and  Pressure 

The  airblast  energy  transferred  to  the  surrounding  medium  from  vessel  rupture 
depends  on  the  internal  temperature  and  pressure  of  the  gas  at  burst.  Fill  gas 
temperatures  were  generally  about  330°K,  although  for  test  7,  they  climbed  to 
37^°K.  At  these  low  temperatures,  vessel  pressure  at  rupture  did  not  significantly 
affect  airblast,  i.e.,  airblast  pressures  from  the  50,000-psi  ruptures  were  not 
measurably  higher  than  those  from  the  15,000-psi  ruptures.  The  reason  for  this  is 
that  there  is  very  little  difference  in  the  expansion  energy  of  argon  at  these 
pressures.  Consider  the  expansion  data  in  Table  3.1  below  based  on  a real  argon 
equation  of  state  from  reference  (U).  The  data  are  for  .a  1 ft3  volume  of  argon 

(!+)  Din,  F.  "Thermodynamic  Functions  of  Gases,  Volume  2:  Air,  Acetylene,  Ethylene, 

Propane,  and  Argon",  Butterworths,  London,  1956 
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at  290°K.  The  expansion  energy,  i.«. the  work  done  against  the  atmosphere,  is 
taken  as  the  difference  in  internal  energy  after  isentropic  expansion  to  a pressure 
of  1 atmosphere. 

TABLE  3.1  ENERGY  IN  ARGON  PRESSURE  VESSELS 


Vessel 
Pressure 
in  psi 

Density 
in  gm/cc 

Internal  Energy 
drop  on  expansion 
cal/gm 

Total  Energy  for 
1 Ft ^ Volume 
cal 

TNT  Equivalent 
in  pounds  for 
1 Ft3 

15,000 

0.99 

13.1*8 

3.81  x 105 

0.82 

30,000 

1. 2U 

12.1*5 

1*.38  x 105 

0.95 

50,000 

1.39 

11.81 

!*.61*  x 105 

1.00 

* Based  on  heat  of  detonation  of  1018  cal/gram 


There  is  little  difference  in  the  expansion  energies  at  these  pressure  levels.  The 
shell  of  a 50,000-psi  vessel  is  about  1+  times  as  massive  as  that  of  a 15,000-psi 
vessel  and  therefore  absorbs  relatively  more  kinetic  energy  than  does  the  lighter 
vessel.  This  should  decrease  the  differences  even  further. 

The  remarks  on  energy  equivalence  in  the  remainder  of  this  discussion  do  not 
imply  that  tank  rupture  airblast  is  identical  to  blast  generated  by  equivalent  TNT 
charges.  The  TNT  equivalence  given  in  Table  3.1  is  used  as  an  energy  unit.  As 
may  be  seen  in  Section  3.5,  neither  the  overpressure  versus  distance  nor  the  impulse 
versus  distance  curves  for  tank  rupture  are  identical  to  those  from  TNT.  While  the 
energy  available  from  vessel  rupture  does  not  increase  above  that  for  an  ideal  gas, 
shock  wave  parameters  are  still  a strong  function  of  gas  temperature.  Figure  E.  6 
shows  the  dependence  of  the  initial  air  shock  pressure  on  temperature  for  real  and 
ideal  argon  for  several  rupture  pressures. 

Consider  now  the  case  where  the  vessel's  internal  gas  temperature  is  very  high. 
As  temperature  is  increased  (somewhere  above  1350°K  for  15,000-psi  pressure)  argon 
behavior  approaches  that  of  an  ideal  gas.  Since  for  an  ideal  gas,  the  expansion 
energy  is  independent  of  temperature  and  molecular  weight,  the  ideal  gas  internal 
energy  represents  an  upper  limit  to  the  blast  energy  available  from  an  argon  vessel 
rupture. 

This  upper  limit  is  shown  in  Figure  3.8.  Note  that  for  1 ft 3 of  arson  burst  at 
15,000  pri.  and  1273°K,  the  TNT  energy  equivalent  approaches  2 pounds.  Further  fill  sas 
temperature  increases  will  not  increase  this  energy  beyond  2.15  lbs  TNT.  Argon 
approaches  an  ideal  monotomic  gas  asymptotically  with  temperature  increase.  The 
reason  for  this  is  that  at  constant  pressure,  a temperature  increase  is  accompanied 
by  a density  decrease.  At  low  densities,  the  intermolecular  forces  that  cause 
non-ideality  become  unimportant.  At  pressures  above  15,000-psi,  temperatures 
greater  than  1300°K  are  required  for  argon  to  approach  ideal  behavior.  Isentropes 
generated  from  the  data  in  refei’ence  (5)  are  shown  in  Figure  3.9.  From  these, 
argon  must  be  near  2100°K  at  50,000-psi  to  approach  ideal  behavior. 

(5)  Brahinsky,  H.S.  and  Neel,  C.A.,  "Tables  of  Equilibrium  Thermodynamic  Properties 
of  Argon  Vol.  Ill  Constant  Entropy".  Arnold  Engineering  and  Development  Co.,  No. 
AEDC-TR  69-19,  Vol.  Ill  AD6814532,  March  1969. 
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The  upper  limits  for  other  gases  are  shown  in  Figure  3.10.  The  upper  energy 
bound  for  a diatomic  gas  at  15*000  psi  is  about  3.25  lbs  TNT  equivalent  versus 
about  2.15  lbs  equivalent  for  argon. 

The  free  air  overpressures  calculated  via  WUNDY,  a one-dimensional  hydrocode 
and  plotted  in  Figure  3.11  illustrate  the  dependence  of  airblast  pressure  on  argon 
temperature. 

3.3  Effects  of  Volume 

Of  all  the  factors  that  determine  airblast  from  pressure  vessel  rupture,  the 
effect  of  change  in  vessel  size  may  be  easiest  to  handle.  Cube  root  scaling,  as 
used  to  compare  airblast  parameters  in  explosives  studies,  seems  applicable  here. 

Cube  root  scaling  is  most  applicable  to  spherical  or  point  source  explosions 
in  free  air.  Change  the  charge  shape,  put  a case  around  the  charge,  or  fire  the 
charge  near  a boundary  and  the  use  of  cube  root  scaling  must  be  qualified. 


Insufficient  vessel-burst  data  exists  to  verify  the  use  of  cube  root  scaling. 
Certainly,  one  expects  it  to  apply  to  spherical  tanks  burst  under  the  conditions 
where  the  only  variable  is  the  vessel  volume.  Thus  distances,  times,  and  shockwave 
impulses  that  characterize  the  airblast  from  vessels  of  different  sizes  are 
related  by  the  ratio  of  the  cube  root  of  their  volume. 

There  are  more  uncertainties  for  vessels  not  spherical  in  shape.  The  most 
likely  vessel  shape  will  be  cylindrical.  Conventional  explosive  cylinders  generate 
an  airblast  field  near  the  charge  that  reflects  the  charge  dimensions  and  the 
detonation  scheme  (reference  6).  However  in  the  case  of  a cylindrical  charge, 
asymmetries  in  the  blast  front  disappear  by  the  time  the  blast  wave  has  propagated 
to  the  point  where  the  airblast  overpressures  are  about  5 psi.  From  then  on,  the 
front  expands  spherically. 

Perturbations  due  to  the  detonation  scheme  of  a conventional  explosion  never 
appear  on  the  tank  rupture  blast  front.  These  are  replaced  by  perturbations  caused 
by  the  vessel  failure  scheme.  The  high  initial  shock  velocities  that  allow  the 
shock  front  from  a conventional  explosion  to  approach  spherical  symmetry  at  late 
times  may  not  exist  in  the  vessel  rupture  situation.  Therefore,  the  blast  front 
from  vessel  rupture  may  never  overcome  its  initial  asymmetries.  In  spite  of  these 
problems,  no  better  method  exists  to  relate  the  blast  from  different  sized  vessels 
than  one  using  cube  root  scaling.  As  the  vessel  rupture  approaches  an  ideal 
spherical  explosion,  cube  root  scaling  will  obtain. 


3.h  Effects  of  a Confining  Structure 

Most  working  high  pressure  vessels  will  be  housed  in  some  kind  of  structure 
built  to  contain  or  modify  the  blast  and  fragment  patterns  should  accidental 
rupture  occur.  Estimates  of  damage  to  personnel  and  nearby  structures  must  take 
into  account  the  effect  of  this  confining  structure  on  damaging  airblast  pressures. 

1 

(6)  Wisotski,  J.  and  Snyer,  H.  "Characteristics  of  Blast  Waves  Obtained  from  Cylin- 
drical High  Explosive  Charges".  Denver  Research  Institute,  DR1^2286,  Nov  1965 
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The  confining  structure  test  described  in  Annexes  A and  C show  some  of  the 
effects  resulting  from  the  presence  of  a five-sided  enclosure  about  a bursting 
high  pressure  vessel.  Curves  fitted  to  the  measured  data  from  Figure  C.17  are 
shown  as  solid  lines  in  Figure  3.12. 

For  the  confined  test  on  which  the  solid-line  curves  in  Figure  3.12  are  based, 
the  vessel  burst  at  a point  that  directed  the  major  gas  motion  against  a sidewall 
of  the  structure. 

If  the  vessel  had  burst  in  a manner  so  as  to  direct  the  argon  jet  out  through 
the  open  wall  of  the  structure,  then  airblast  overpressures  would  have  been  even 
higher.  The  dashed  line  in  Figure  3.12  was  obtained  from  the  average  of  the  curves 
in  Figure  3.5.  This  curve  probably  represents  the  maximum  airblast  overpressures 
in  front  of  the  structure. 

The  dotted  curve  in  Figure  3.12  represents  the  average  of  the  curves  shown  in 
Figure  3.6.  This  curve  gives  pressures  along  lines  up  to  71°  from  the  argon  Jet. 

Even  for  these  conditions,  airblast  overpressures  at  some  distance  from  the  burst- 
ing vessel  are  higher  than  would  have  been  predicted  for  a confined  burst.  Thus, 
under  certain  burst  conditions  the  five-sided  confining  structure  does  not 
reduce  airblast  overpressure. 

Pressure  on  the  interior  walls  will  also  depend  on  burst  geometry  and 
internal  gas  temperature.  The  2-dimensional  Hydrocode,  TUTTI,  will  calculate  the 
entire  pressure  load  on  the  walls  as  a function  of  time  for  all  burst  conditions. 

3.5  Comparison  of  Data 

Measured  airblast  overpressures  are  compared  vith  calculated  data  in  Figure 
3.13.  The  upper  bound  of  the  shaded  area  was  arrived  at  by  taking  the  average  of 
the  data  from  Figure  3.5;  the  lower  bound  is  the  average  of  the  data  from  Figure 
3.7. 

The  data  from  the  one-dimensional  calculation,  WUNDY , was  modified  to  fit  the 
burst  geometry  of  these  experiments.  Specifically,  data  above  1*0  psi  are  from 
Figure  E.ll.  Data  below  20  psi  was  obtained  from  Figure  3.11.  In  this  case,  the 
WUNDY  Curve  for  all  burst  pressures  with  an  argon  temperature  of  290°K  (17°C)  was 
used.  The  one-dimensional  calculations  are  for  a free  air  burst.  To  account  for 
the  presence  of  the  ground,  the  energy  output  was  doubled.  Therefore,  distances 
were  multiplied  by  (2)  or  1.26.  The  TNT  curve  was  obtained  from  Figure  C.l 
using  the  equivalent  weight  data  from  reference  7. 

Airblast  overpressures  at  a ground  range  of  1.0,  Figure  3.13,  range  from  10 
psi  for  the  lower  bound  of  the  measured  data  to  about  1,000  psi  for  the  TNT  data. 

The  curves  converge  at  lower  overpressure  levels.  The  upper  bound  of  the  measured  data 


(7)  Swisdak,  M. , ''Explosion  Effects  and  Properties:  Part  1 - Explosion  Effects  in 

Air",  Naval  Surface  Weapons  Center  Report  NSWC/WOL/TR  75-116 
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from  the  TUTTI  claculation,  and  the  TNT  data  coincide  at  overpressures  below  about 
5.0  psi.  Thus  either  the  TUTTI  results  or  the  TNT  curve  represent  an  upper  bound 
to  the  airblast  overpressures  generated  by  vessels  ruptured  in  this  investigation. 

1*.  FRAGMENT  DATA 


Detailed  fragment  data  are  given  in  Annex  D.  All  vessels  burst  into  two  pieces. 
Of  these  lU  pieces,  10  were  recovered.  Weights  of  the  recovered  sections 
ranged  from  51.1  pounds  to  271  pounds.  The  recovered  sections  were  those  that  hit 
the  arena  walls  and  were  stopped  or  those  that  hit  the  ground  or  firing  pad  and 
traveled  only  a short  distance.  The  positions  where  the  pieces  of  the  vessels 
were  found  are  shown  in  Figure  D.l. 

The  fragment  velocity  system  did  not  function  satisfactorily.  Either  the 
pieces  of  the  vessel  failed  to  hit  a screen;  or  when  they  did,  the  electronic 
counters  did  not  function.  Fortunately,  a piece  of  the  vessel  hit  the  firing  pad 
on  all  the  unconfined  tests  except  shot  5.  This  produced  a distinctive  signal 
on  the  time  resolved  pressure  gage  recordings.  The  points  where  the  vessel  seg- 
ments hit  the  firing  pad  and  the  original  position  of  the  vessel  were  used  to 
determine  the  distance  the  fragment  traveled.  Early  fragment  velocities  were 
determined  from  this  information.  For  the  test  where  the  vessel  half  did  not  hit 
the  firing  pad,  shot  5,  (see  Table  D.l)  one  section  hit  velocity  screens. 

The  mean  velocity  measured  for  the  approximately  50-pound  portion  of  the 
15,000-psi  vessels  was  320  feet/second.  The  minimum  velocity  measured  was  310  feet/ 
second  and  the  maximum  was  330  feet/second. 

For  the  30,000-psi  vessel,  the  mean  velocity  measured  was  353  feet/second  for 
the  1U7  pound  segment  of  shot  3 and  27^  feet/second  for  the  108-pound  segment  of 
shot  7.  The  average  of  all  fragment  velocities  measured  on  shots  3 and  7 is  326 
feet/second  with  the  highest  value  measure  373  feet/second  and  the  minimum  velocity 
measured  250  feet/second.  The  wide  velocity  differences  are  between  the  two  shots 
rather  than  between  individual  measurements  on  a single  shot. 

For  the  50,000-psi  vessels,  the  average  velocity  measured  for  thv.  258  pound 
portion  of  the  vessel  used  on  shot  ^ was  215  feet/second.  The  only  velocity  screen 
measurements  were  made  on  shot  5.  The  heavier  portion  of  the  vessel,  271  pounds, 
tiaveled  17*5  feet  at  an  average  velocity  of  273  feet/second.  The  lighter  portion, 
ll*2  pounds  traveled  a distance  of  6.8  feet  in  the  opposite  direction  at  an  average 
velocity  of  278/feet/second.  Average  velocity  for  all  50,000-psi  Vessel  fragments 
was  235  feet/second. 

Case  velocities  from  the  one-dimensional  code  calculation  are  shown  in  Figure 
U.l  for  the  15,000-psi  vessel.  Note  that  measured  velocities  are  less  than  half 
the  calculated  value.  A method  for  calculating  case  velocity  from  reference  3 gives 
similar  results,  i.e.,  higher  than  measured  values  by  a factor  of  2.  No  explanation 
for  this  discrepancy  is  readily  available.  A sophisticated  technique  for  fragment 
velocity  calculation  is  given  in  reference  8.  However  it  was  not  applied  to  these 
experiments. 


Baker,  W.  et.  al.  Assembly  and  Analysis  of  Fragmentation  Data  for  Liquid 
Propellant  Vessels,"  NASA  CR  13U538  Prepared  for  NASA  Lewis  Research  Center, 

Aerospace  Safety  Research  and  Data  Institute  by  Southwest  Research  Institute,  San 
Antonia,  Texas,  Jan  1971* 

UNCLASSIFIED 


UNCLASSIFIED 

NSWC/WOL/TR  75-8? 


5.  SUMMARY 

Several  factors  control  the  airblast  field  from  the  pneumatic  hurst  of  high 
pressure  vessels.  Those  related  to  the  vessel  and  its  contents  include  the  type, 
temperature,  volume,  and  pressure  of  the  gas  in  the  vessel  and  the  geometry  of 
rupture. 

For  a vessel  filled  with  a monotomic  gas  and  burst  at  a temperature  below 
about  330° K,  increasing  burst  pressure  from  15,000  to  50,000  psi  does  not  measurably 
increase  blast  energy.  However,  at  elevated  fill  gas  temperatures,  rupture 
energy  is  increased,  'ihe  upper  limit  is  reached  when  internal  temperatures  are 
high  enough  to  cause  the  fill  gas  to  approach  ideal  gas  behavior.  These  temperatures 
are  about  1300°K  for  15,000  psi  and  about  2100°K  for  50,000  psi.  Further  increases 
in  fill  gas  temperatures  will  not  significantly  increase  expansion  energy.  Thus  the 
ideal  gas  assumption  sets  an  upper  bound  to  the  expansion  energy  from  high  pressure 
vessel  rupture.  However,  the  way  in  which  this  expansion  energy  increase  influences 
airblast  overpressures  is  non-linear.  Airblast  overpressure-distance  curves  must 
be  computed  for  each  burst  temperature/pressure  regime. 

Blast  output  from  spherical  high  pressure  vessels  of  different  volumes  burst 
at  the  same  temperature  and  pressure  can  be  related  by  dividing  distances  and 
times  by  the  cube  root  of  the  volume.  Airblast  from  cylindrical  vessels  may  be 
handled  similarly  by  considering  them  to  be  an  equivalent  volume  sphere  if  the 
length  to  diameter  ratio  of  the  cylinder  is  not  much  greater  than  2/l.  It  is 
believed  that  this  procedure  is  relatively  accurate;  however,  the  degree  of 
certainty  cannot  be  determined  given  the  paucity  of  available  data. 

The  blast  field  from  pressure  vessel  rupture  is  strongly  influenced  by  burst 
geometry.  The  vessels  used  in  this  investigation  burst  into  two  pieces  creating 
strong  argon  motion  along  a line  defined  by  the  vessels  center  and  the  point  on 
the  shell  where  first  rupture  occurs.  Airblast  overpressures  at  gages  along  the 
line  of  the  Jet  show  overpressures  a factor  of  b or  more  greater  than  pressures 
along  a line  in  the  opposite  direction  from  the  Jet. 

Two  methods  were  used  to  calculate  the  airblast  from  bursting  high  pressure 
vessels.  A two-dimensional  hydrocode,  TUTTI,  was  used  to  calculate  airblast 
overpressures  near  the  vessel  and  in  non-ideal  situations,  i.e.,  in  the  presence 
of  obstructions.  A one-dimensional  hydrocode,  WUNDY,  may  be  used  to  arrive  at 
airblast  overpressures  in  uncomplicated  situations  or  at  distances  greater  than  5 
or  6 vessel  radii  from  the  burst.  TUTTI  tends  to  give  airblast  overpressures  in 
the  high  range  of  the  measurements  — an  acceptable  situation  where  the  data  is  to 
be  used  to  establish  safety  criteria.  WUNDY  tends  to  give  pressures  in  the  middle 
range  of  the  measured  data. 

The  wall  fragments  observed  in  this  test  program  weighed  from  50  pounds  for 
the  15,000  psi  vessel  to  271  pounds  for  the  50,000  psi  vessels.  These  were 
accelerated  to  velocities  up  to  350  feet/second.  For  many  situations,  fragments 
are  a greater  threat  than  is  the  blast  generated  by  high  pressure  vessel  rupture. 
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FIG.  2.3  50,000-PSI  TEST  VESSEL 


21 

UNCLASSIFEID 


GAGE  LINE  1 


FIG.  3.1  ESTIMATED  DIRECTION  OF  ARGON  JETTING 
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ANNEX  A EXPERIMENTAL  PLAN 


A.l  Pressure  Vessels 

All  test  vessels  were  spherical  containers  with  an  internal  radius  of  7-5 
inches  and  a volume  of  one  cubic  foot.  They  were  made  from  preformed  hemispheres 
of  T-l  steel  welded  together.  A section  of  high  pressure  tubing  welded  into 
the  vessel's  boss  provided  a gas  inlet.  This  arrangement  is  shown  in  Figure  2.3. 
Three  of  the  vessels  were  designed  to  burst  at  15,000  psi,  two  were  designed  to 
burst  at  30,000  psi.  and  the  remaining  two  were  designed  to  burst  at  50,000  psi. 

The  yield  strength  of  T-l  steel  is  137  Kpsi,  the  ultimate  strength  is  151  Kpsi. 
Measured  yiild  and  ultimate  strengths  of  the  equatorial  weld  material  were  lib 
Kpsi  and  121  Kpsi  respectively.  Other  vessel  data  pertinent  to  this  investigation 
are  given  in  Table  2.1. 

3 

The  argon  fill  gas  weight  in  Table  2.1  is  given  for  a 1 ft  volume.  It  is 
based  on  the  compression  data  for  argon  given  in  Section  2.2  of  this  report.  A 
standard  cubic  foot  (SCF)  of  argon  as  used  here  is  given  as  0.111  pounds  of 
argon  at  290°K  and  lU . U psiA,  (reference  A.l). 

A. 1.1  Pressure  Vessel  Expansion 

S' 

l 

Pressure  vessel  expansion  was  recorded  as  a function  of  pressure  from 
beginning  of  pressurization  to  burst.  Measurements  were  made  along  two  mutually 
perpendicular  circumferences  as  shown  in  Figure  A.l.  The  system  consisted  of 
cords  wrapped  around  the  vessel  and  attached  to  linear  potentiometers.  An  increase 
in  vessel  circumference  moved  the  slide  on  the  potentiometer  producing  a resis- 
tance increase  proportional  to  the  increase  in  circumference.  Potentiometer  out- 
puts were  recorded  on  chart  recorders. 


A. 2 Pressurization  System 


The  multi-stage  pressurization  system,  shown  schematically  in  Figure  A. 2,  was 
manufactured  by  the  Haskel  Engineering  and  Supply  Company  of  Burbank,  California*. 
The  argon  was  supplied  from  bottles  at  a pressure  of  between  1,000  and  2,350  psig. 
This  pressure  was  increased  to  35,000  psig  by  the  model  lU  ATS-315C  gas  booster. 
The  model  15939-2  intensifier  increased  the  pressure  from  35,000  psig  to  a maximum 
of  60,000  psig.  The  intensifier  piston  was  driven  by  a model  GS-100-C  hydraulic 
pump.  Both  the  booster  and  the  hydraulic  pump  were  driven  by  an  air  compressor 
supplying  200  SCF/minute  at  100  psig. 

Pumping  rate  of  the  system  is  a function  of  the  air  drive  pressure,  the  argon 
supply  pressure,  and  the  outlet  pressure.  Pressurization  times  for  the  vessels 
depended  on  pumping  rate  and  the  pressure  in  the  vessel.  Compression  data  for 
argon  is  shown  in  Figure  A.  3 (reference  A.l).  Average  pumping  rates  were  about 
9 SCF/minute  for  the  15,000_psi  vessels,  5 SCF/minute  for  the  30,000-psi  vessels, 

(A.l)  Din,  F.  fTThermodynamic  Functions  of  Gases,  Volume  2:  Air,  Acetylene, 

Ethylene,  Propane,  and  Argon  "Butterworths,  London,  1956. 

•Mention  of  a manufacturer's  products  by  name  constitutes  neither  criticism  nor 
endorsement  by  the  Naval  Surface  Weapons  Center. 
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and  U . 9 SCF/minute  for  the  50,000-psi  vessels. 

Internal  vessel  pressures  were  measured  with  a transducer  near  the  inlet  of  the 
test  vessel.  As  shown  in  Figure  A. 2 the  gage  was  fed  through  a 20  foot  section 
of  high  pressure  tubing  connected  into  the  main  line  280  feet  from  the  pump  and 
20  feet  from  the  test  vessel.  The  static  error  band  for  the  transducer  was+0.1*# 
full  scale  based  on  a terminal  straight  line.  Thermal  effects  gave  a zero  voltage 
shift  of  0.05#  of  full  scale  per  °C.  The  gage  output  was  recorded  in  real  t^me 
on  a chart  recorder. 

A thermocouple  extending  into  the  test  vessel  continuously  monitored  the  argon 
temperature  in  the  vessel.  The  output  of  the  thermocouple  was  recorded  on  a strip 
chart  recorder.  System  circuitry  and  a typical  calibration  curve  for  a chromel/ 
alumel  device  are  shown  in  Figure  A.l*. 

A. 3 Airblast  Instrumentation 

A. 3.1  Blast  Gages  and  Layout 


Airblast  generated  by  vessel  rupture  was  measured  along  two  lines  for  the  six 
unconfined  tests  as  shown  in  Figure  A. 5.  The  seven  gages  used  on  each  line  were 
deployed  in  the  following  manner.  The  first  two  gcges,  1-1  and  2-1  from  Figure 
A. 5,  were  flush  mounted  in  the  concrete  firing  pad  at  a ground  range  of  1 foot. 
These  gages  were  either  Kulite  Semiconductor  Product's  model  HSK-375  or  Susquehanna 
Instrument's  model  ST-1*.  The  second  two  gages  were  located  at  mid  vessel  height 
and  1 foot  from  the  vessel's  perimeter.  These  were  Model  ST-7  gages  made  by 
the  Susquehanna  Instrument  Company.  The  third  two  gages,  positions  1-3  and  2-3, 
also  mounted  in  the  firing  pad,  were  at  a ground  range  of  2.5  feet.  These  gages 
were  models  LC-70  or  LC-71  gages  made  by  Celesco  Industries.  Positions  1-1*  and 
2-1*,  1-5  and  2-5,  and  1-6  and  2-6  were  located  in  individual  2x2  feet  concrete 
pads  at  ground  ranges  of  5,  10,  and  16.5  feet  respectively.  The  gages  used  in 
these  positions  were  also  model  LC-70's  or  LC-71's.  Positions  1-7  and  2-7  were  at  a 
ground  range  of  60  feet.  These  gages,  model  LC-33's  made  by  Celsco  Industries,  were 

mounted  1 foot  above  the  ground.  Another  LC-33  gage,  position  3-1  in  Figure  A. 5, 
was  mounted  1 foot  above  the  ground  and  at  a ground  range  of  13.2  feet. 

Figure  A. 6 shows  the  gage  layout  for  the  confined  test.  Changes  from  the  un- 
confined test  layout  were:  the  closest  gage  positions  were  removed,  an  additional 

gage  line  was  installed,  and  six  face-on  gages  were  installed  flush  with  the  walls 
of  the  confining  structure.  These  gages  were  supported  by  mounts  outside  the 
structure,  as  shown  in  Figure  A. 7.  A potting  compound  isolated  the  gage  mounts 
from  the  structure  walls  to  prevent  unwanted  acceleration  being  imparted  to  the 
gages  by  the  wall  motion.  (Figure  A.8). 

One  ST-1*  and  one  LC-70  gage  was  used  in  each  wall.  A thin  sheet  of  pyrolytic 
graphite  was  glued  to  the  face  of  each  LC-70  gage  to  insulate  the  gage  against  the 
low  temperature  of  the  expanding  argon. 

A.  3.2  Blast  Instrumentation 


Gage  signals  were  fed  to  charge  amplifiers  located  at  the  gages.  Signals  from 
the  charge  amplifiers  were  fed  through  coaxial  cables  to  the  recording  equipment 
located  in  an  instrument  trailer  some  500  feet  from  ground  zero. 
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The  recording  system  consisted  of  12  oscilloscope  channels  and  lL  magnetic 
tape  channels.  The  oscilloscopes,  tape  recorder,  and  the  auxiliary  equipment  used 
for  time  bases,  synchronization,  and  power  were  housed  in  a special  instrumentation 
trailer.  This  and  similar  systems  have  been  used  at  the  Center  over  the  past 
several  years.  They  are  described  in  reference  A. 2. 

Two  of  the  tape  recorder's  channels  were  used  to  record  time  code  information 
and  playback  signals.  Thus  12  tape  channels  and  12  oscilloscope  channels  were 
available  for  blast  recording.  Since  15  gages  were  used  for  the  unconfined  tests, 

9 gages  were  recorded  on  both  tape  and  oscilloscopes.  For  the  confined  test,  21 
gages  were  used  so  that  only  three  channels  were  recorded  on  both  systems. 

A.  1+  Confining  Structure 

One  15,000  psi  vessel  was  burst  inside  a 5-sided, cubical  steel  structure. 
Figures  A. 7 and  A. 8.  The  open  side  of  the  vessel  was  perpendicular  to  the  ground. 
The  test  vessel  was  centrally  located  in  the  structure. 

The  structure's  internal  volume  was  determined  by  roughly  scaling  it  to  the 
volume  of  an  autoclave  facility  in  operation  at  the  Oak  Ridge  Y-12  Plant.  The 
autoclave  facility's  volume  is  about  lit, 700  ft^.  The  autoclave  external  volume 
is  about  500  ft"5  giving  a volume  ratio  of  29.  The  test  vessel's  volume  is  1.2 
ft 3.  Thus  the  confining  structure's  volume  of  35  ft^  roughly  scaled  the  Y-12 
Plant  Facility. 

A. 5 Fragment  Recovery  and  Velocity  System 

A.  5.1  Fragment  Recovery 

The  fragment  recovery  and  velocity  system  design  was  based  on  the  possibility 
that  the  test  vessels  might  rupture  in  a brittle  manner,  i.e.,  break  up  into 
several  fragments. 

A semi-circular  arena  was  constructed  to  recover  a sample  of  the  vessel 
fragments.  The  arena  was  made  from  Celotex  bundles  placed  as  shown  in  Figure  2.2. 
Each  bundle  was  made  from  81+  sheets,  each  1+'  x 8'  x 0.5"  to  form  a bundle  1+'  x 
8'  x 1+2".  These  bundles  were  deployed  as  shown  in  Figure  2.2  to  form  the  arena. 

A.  5.2  Fragment  Velocity  System 

All  fragment  velocity  measurements  were  made  using  arrival  time  techniques. 

The  instrumentation  consisted  of  seven  electronic  counters.  All  counters  were 
started  by  a zero  time  signal.  This  signal  was  generated  by  the  interruption  of 
an  electrical  current  in  a break  wire  wrapped  around  the  vessel's  outer  surface. 

The  counters  were  stopped  by  "Velocity  screens"  placed  at  various  positions  about 
the  vessel.  These  screens  consisted  of  a continuous  wire,  wrapped  about  a frame- 
work. The  wire  carried  a current  that  was  interrupted,  when  severed  by  a fragment, 
to  supply  a stop  signal.  Signal  shapers  were  used  between  the  break  wires  and 
the  counters  to  insure  uniform  signals  to  the  counter  inputs. 

A. 2 Tussing,  R.,  "A  Four  Channel  Oscilloscope  Recording  System",  Naval  Ordnance 
Laboratory,  N0LTR  65-21,  18  May  1965 


A_3 

UNCLASSIFIED 


UNCLASSIFIED 

NSWC/WOL/TR  75-87 


Screen  locations  are  shown  in  Figures  2.1  and  2.2.  Screen  sizes  varied 
according  to  their  location  as  may  be  seen.  The  entire  arena  wall  was  divided  up 
into  screens  1 and  2.  Screens  3,  1»,  and  5 were  located  on  a teepee  arrangement 
above  the  vessel.  Screen  6 was  located  behind  the  vessel  and  Screen  7 was  located 
around  the  test  vessel.  Because  the  vessel  failed  to  fragment  in  a brittle  manner, 
the  screen  arrangement  was  varied  for  the  seven  tests. 
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ANNEX  B VESSEL  RUPTURE  HISTORIES 


This  Section  is  included  to  display  the  internal  vessel  pressure  and  tempera- 
ture, and  the  vessel  expansion  history  recorded  during  pressurization.  Figure  B.l 
shows  a 15,000-psi  vessel  on  the  stand  ready  for  test.  In  this  case  the  equational 
vessel  veld  plane  is  parallel  with  the  ground/ firing  pad.  Since  the  rupture 
geometry  greatly  affected  the  airblast  field,  a description  of  each  vessel  rupture 
is  given. 

B.l  Rupture  Histories  of  15.000-PBi  Vessels 

Internal  vessel  pressures  and  temperatures  are  shown  for  the  three  15,000-psi 
bursts  in  Figures  B.2,  B.3,  and  B.l*.  No  preburst  vessel  expansion  was  indicated 
on  any  15,000-psi  vessel. 

Only  actual  pumping  times  are  plotted  in  the  Figures.  In  some  cases,  ergon 
was  bled  into  the  test  vessel  directly  from  the  supply  bottles.  This  bleed-time 
was  not  included  in  the  pressurization  times  given  for  any  vessel.  A total 
pumping  time  of  1+5  minutes  was  required  to  pressurize  the  shot-1  vessel  from 
2150  psi  to  burst  at  lU,950  psi  as  shown  in  Figure  B.2. 

The  pump  was  halted  several  times  during  pressurization  to  perform  various 
tasks  such  as  changing  argon  supply  bottles,  clearing  the  firing  area,  or 
adjusting  instrumentation.  These  halts  are  marked  by  discontinuous  changes  in  the 
pressurization  or  temperature  curves.  After  a pause,  the  resumption  of  pumping 
was  accompanied  by  a sharp  temperature  increase.  Note  those  in  Figure  B.2  at 
0,  12,  and  25  minutes.  Beyond  30  minutes,  the  system  was  slowly  approaching  equil- 
ibrium temperature.  Net  temperature  rise  for  shot  1,  Figure  B.2  was  35°C. 

Figure  B.3  shows  the  pressurization  rate  end  internal  gas  temperature  for 
shot  2.  The  pumping  began  when  the  internal  vessel  pressure  was  500  psi  and 
continued  until  the  rupture  occurred  at  ll+,750  psi  after  some  76  minutes  pumping 
time.  The  pumping  rate  was  very  slow  up  to  5,000  psi  . Therefore,  the  internal 
gas  temperature  decayed  following  the  initial  Jump  at  the  resumption  of  pumping 
after  a pause.  The  ambient  temperature  remained  constant  during  the  pumping  phase; 
the  internal  argon  temperature  increased  from  5 to  27.5°C. 


Figure  B.l*  shows  the  pressurization  rate  and  internal  gas  temperature  for 
shot  6,  the  confined  vessel  test.  This  vessel  ruptured  at  ll+,l+00  psi  after  73 
minutes  pumping.  Internal  gas  temperature  rose  from  21°C  ambient  to  65°C,  a 
net  increase  of  1+1+°C. 

Figure  B.5  shows  the  estimated  rupture  geometry  for  the  shot-1  vessel.  Prior 
to  rupture,  the  vessel  was  set  up  as  shown  in  Figure  B.l,  in  this  case  with  the 
equatorial  weld  plane  parallel  with  the  earth/firing  pad  surface.  At  rupture, 
the  vessel  opened  like  a clam  shell  with  the  first  opening  occuring  along  the 
equatorial  weld  at  a point  intersected  by  the  arrow  shovinc  the  shot-1  Jet  direc- 
tion. After  opening,  the  reaction  forced  the  bottom  half  downward  and  away  from 
the  Jet  direction.  It  then  hit  the  pad  and  ricocheted  upward  over  the  arena  walls 
and  landed  some  665  feet  from  ground  zero.  Based  on  the  trajectory  of  the  bottom 
half,  the  shot-1  vessel's  top  half  rocketed  off  at  an  elevation  angle  of  5^*°  in  the  same 
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direction  as  the  bottom  half.  (See  Figure  D.l). 

A similar  scenario  seems  to  fit  the  rupture  geometry  of  the  shot-2  vessel. 

It  too  burst  at  the  equator  but  with  rupture  beginning  at  a point  in  the  direction 
of  the  estimated  argon  Jet  for  shot-2  shown  in  Figure  B.5.  The  bottom  half  was 
forced  downward,  striking  the  firing  pad  and  ricocheting  upward  and  outward.  The 
bottom  half  was  found  552  feet  from  ground  zero.  Based  on  the  trajectory  of 
the  bottom  half,  the  top  half's  trajectory  was  estimated  to  have  an  elevation 
angle  of  6L°  and  to  be  directly  away  from  the  estimated  shot-2  argon  Jet  (Figures 
B.5  and  D.l);  the  top  half  was  not  found. 

The  remaining  15,000-psi  vessel  was  burst  in  the  confining  structure.  In  this 
case,  the  equatorial  weld  plane  was  perpendicular  to  the  earth/firing  pad  surface. 
The  bottom  half  (the  half  containing  the  boss  and  gas  inlet  tube)  faced  toward 
the  opening  of  the  confining  structure  as  shown  in  Figure  A. 7. 

This  vessel  also  ruptured  at  the  equatorial  weld.  The  top  half  of  the  vessel 
struck  the  rear  wall  of  the  confining  structure.  The  impact  mark  was  centered 
3"  below  and  7"  to  the  left  of  the  rear  wall's  center.  The  bottom  or  outward 
facing  half  of  the  vessel  struck  the  ground  along  a trajectory  that  was  the  mirrow 
image  of  the  top  half's  path. 

This  opening  pattern  directed  an  argon  Jet  toward  a point  on  the  right  wall 
(facing  the  opening  and  looking  in)  of  the  structure.  The  resultant  forces  turned 
the  structure  31*0  clockwise  and  moved  it  back  some  two  feet.  (Figure  B.6). 

B.2  Rupture  Histories  of  30.000-ps1  Vessels 

Rupture  histories  recorded  for  the  30,000-psi  vessels  are  shown  in  Figure 
B.7  and  B.8.  The  shot-3  vessel  burst  at  3^, ^00  psi  after  98  minutes  pumping 
time  (Figure  B.7).  Pumping  began  with  the  vessel  at  2300  psi.  The  internal 
vessel  temperature  rose  from  ll*°C  ambient  to  23°C  within  2 minutes  after  pumping 
began.  It  then  began  a decrease  for  the  next  16  minutes  because  of  the  slow 
pressurization  rate.  At  this  point,  fresh  gas  supply  bottles  increased  the  pumping 
rate  and  the  rate  of  temperature  rise.  The  temperature  eventually  rose  to  1+3°C 
followed  by  a decrease  during  the  last  four  minutes  of  pumping.  This  decrease 
was  caused  by  a decrease  in  pressurization  rate  as  vessel  pressure  approached 
35,000  psi,  the  limit  of  the  pressurization  system's  first  stage. 

An  increase  of  the  vessel's  circumference  was  measured  as  the  pressure  increased 
from  30,000  psi  to  burst.  This  increase  amounted  to  0.10  inches  on  one  gage  and 
O.065  inches  on  the  other. 

Figure  B.8  shows  the  pressurization  history  for  the  shot-7  vessel.  This 
vessel's  walls  were  machined  down  by  0.12  inches  in  an  effort  to  induce  a failure 
mode  not  connected  with  the  equatorial  weld.  After  115  minutes  pumping  time,  the 
vessel  ruptured  at  31,800  psi  . Because  of  air  compressor  problems,  the  pumping 
system's  intensifier  was  needed  to  complete  pressurization  from  30,000  psi  to 
rupture.  This  accounts  for  the  sawtooth  appearance  of  the  pressure  and  temperature 
curves  above  30,000  psi  . The  intensifier  was  manually  cycled  at  a rate  of  about 
0.5  cycles/minute. 
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The  internal  gas  temperature  reached  a maximum  of  91°C  before  cooling  as  the 
pressurization  rate  fell  with  the  use  of  the  intensifier.  Ambient  temperature 
also  increased  during  pressurization  so  that  the  maximum  net  temperature  rise 
was  72°C. 

The  circumferential  expansion  was  0.43"  and  0.45"  as  measured  by  gage3  1 and 
2 respectively. 

The  shot  3 vessel  was  burst  with  the  equatorial  weld  parallel  with  the  earth/ 
firing  pad  surface.  The  vessel  ruptured  along  the  equatorial  weld.  At  rupture, 
the  bottom  half  of  the  vessel  was  driven  almost  straight  down,  ricocheted  off  the 
firing  pad,  went  up  in  the  air  and  struck  the  arena  wall  on  its  way  down.  The  top 
half  of  the  vessel  was  never  found.  Geometry  of  rupture  and  argon  Jetting  is 
shown  in  Figure  B.9. 

The  shot-7  vessel  was  oriented  with  the  equatorial  weld  perpendicular  to  the 
earth/firing  pad  surface.  The  bottom  or  boss/inlet  half  of  the  vessel  pointed 
away  from  the  arena  along  a line  bisecting  the  angle  between  gage  lines  1 and  2. 

The  shot-7  vessel  also  burst  at  the  equatorial  weld  even  though  the  vessel's 
shell  had  been  thinned  by  0.12  inches.  (The  weld  area  was  not  cut  down).  Both 
the  top  and  bottom  halves  were  found;  the  top  half  after  striking  the  arena  wall 
and  remaining  in  the  arena,  the  bottom  half  after  plowing  a furrow  in  the  dirt 
and  landing  180  feet  from  ground  zero.  Shot-7  rupture  geometry  is  shown  in 

Figure  B.9. 

B.  3 Rupture  Histories  of  50,000-psi  Vessels 

The  parameters  measured  during  pressurization  of  the  50,000-psi  vessels  are 
shown  in  Figure  B.10  for  shot  1*  and  Figure  B.ll  for  shot  5.  Both  vessels  ruptured 
at  50,400  psi  after  about  150  minutes  of  pumping.  The  sawtooth  shape  of  the 
pressure  and  temperature  readouts  in  Figures  B.10  and  B.ll  begin  at  pressures  above 
30,000  psi  when  it  became  necessary  to  use  the  manually  cycled  intensifier.  The 
cycling  rate  is  roughly  0.5  cycles /minute. 

The  temperature  drop  between  cycles  is  seen  to  be  greater  for  the  shot-5 
vessel  in  Figure  B.ll  than  for  the  shot-4  vessel  in  Figure  B.10.  The  net  tempera- 
ture rise  was  36°C  for  shot  4 and  44°C  for  shot  5. 

A circumferential  expansion  of  0.70  inches  was  recorded  on  gage  1 for  shot  4 
(Figure  B.10).  The  other  gage  did  not  register  on  this  vessel.  Gage  2 recorded 
an  expansion  of  0.35"  cm  shot  5.  Gage  1 aid  not  register  for  this  event. 

Both  50,000-psi  vessels  were  placed  on  the  firing  stand  with  the  equatorial 
weld  perpendicular  to  the  earth/firing  pad  surface.  The  bottom  half  of  the  vessel 
pointed  away  from  the  arena  wall  and  along  a line  bisecting  the  angle  between  gage 
lines  1 and  2. 

Both  vessels  were  undercut  so  they  were  thinner  in  the  area  around  the  boss. 
(See  Figure  2.3).  In  both  cases,  failure  started  in  this  area  and  proceeded  in 
a plane  roughly  perpendicular  to  the  equatorial  weld.  Thus  the  vessel  burst  into 
two  pieces  of  unequal  size  as  may  be  seen  in  Figure  D.l. 

B-3 

UNCLASSIFIED 


- 1 jMjfatfcaty  S.faknMiAy-jifiiiirr' 


UNCLASSIFIED 

NSWC/WOL/TR  75-87 

In  the  case  of  shot  1*,  the  rupture  plane  was  nearly  parallel  to  the  ground/ 
firing  pad  surface  with  about  a 10°  tilt  to  the  right  looking  from  the  bottom 
or  boss  half.  The  larger  portion  traveled  downward  and  struck  the  firing  pad. 

It  then  ricocheted,  striking  the  arena  wall  and  remained  within  the  arena. 

(Figure  B.12).  The  upper  portion  (the  smaller  of  the  two)  was  not  found. 

The  she' -5  vessel's  rupture  plane  was  nearly  vertical,  i.e.,  the  plane  tilted 
about  16°  from  the  vertical  in  a counterclockwise  direction  when  viewed  from  the 
boss  end  of  the  vessel.  The  smaller  portion  traveled  upward  and  struck  the  right 
arena  wall  6 feet  above  the  ground.  The  larger  portion  struck  the  ground  3.25 
feet  from  ground  zero.  From  there,  it  bounced  into  the  left  arena  wall  near  the 
ground . 

The  estimated  directions  of  the  argon  Jets,  based  on  the  motions  of  the 
fragments  are  shown  in  Figure  B.12. 
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FIG.  B.9  VESSEL  RUPTURE  GEOMETRY  OF  SHOT  3 AND  SHOT  7;  30,000-PSI. 
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ANNEX  C AIRBLAST  RESULTS 

C.  1 Airblasb  from  Pentollte  Spheres 

The  energy  available  from  the  pressure  vessels  used  in  this  investigation  was 
roughly  equivalent  to  that  available  from  1 pound  of  TNT;  but,  the  airblast 
generation  by  conventional  explosions  and  tank  rupture  arr>  not  comparable.  Since 
the  blast  performance  of  TNT  and  pentolite  are  reasonably  well  known  (see 
reference  C.l  for  instance  ) , it  was  felt  prudent  to  obtain  data 

under  the  conditions  of  this  experiment  to  serve  for  rough  comparisons  and  to  check 
the  blast  instrumentation.  Pentolite  was  used  rather  than  TNT  because  of  the 
difficulty  in  detonating  1-pound  TNT  spheres.  Pentolite  data 
is  easily  scaled  to  TNT  data  using  the  equivalent  weight  information  in  reference 
C.l. 


Two  bjQ-gm  pentolite  spheres  were  detonated  at  a height-of-burst  (HOB)  of 
19-2  inches,  the  nominal  HOB  for  the  test  vessels.  The  same  blast  instrumentation 
in  the  same  arrangement  was  used  for  the  pentolite  tests  and  the  unconfined  vessel 
tests . 

The  pentolite  data  are  tabulated  in  Table  C.l.  Blast  overpressures  are 
compared  with  the  standard  reference  C.l  data  in  Figure  C.l.  vigures  C.2  and  C.3 
show  similar  comparisons  for  positive  shockwave  durations  and  positive  shockwave 
impulses. 

A segment  of  the  standard  free  air  pentolite  curve  is  included  for  comparison 
with  our  one-point  measurements  in  Figure  C.l.  Agreement  between  measured  and 
standard  data  is  satisfactory.  Pressure-time  histories  from  the  pentolite  spheres 
are  shown  in  Figure  C.h. 

C.2  Airblast  from  15 ,000-psi  Vessel  Rupture 

Measured  airblast  data  from  the  two  15,000-psi  vessel  ruptures  are  tabulated 
in  Table  C.2.  Airblast  overpressures  versus  ground  range  have  been  plotted  in 
Figure  C.5. 

The  data  displayed  in  Figure  C.5  is  typical  of  vessel-rupture  data.  At  burst, 
an  asymmetrical  blast  front  develops  as  the  fill  gas  jets  out  at  vessel-wall 
openings.  If  the  vessel  bursts  into  many  pieces,  there  are  many  jets  which  tend 
to  coalesce  into  a spherically  (or  hemispherically ) expanding  shock  front.  The 
pressure  in  the  shock  front  has  a pronounced  effect  on  the  rate  at  which  the  front 
tends  to  coalesce.  This  comes  about  because  at  higher  pressures  the  differences 
in  shock  velocities  between  that  higher  pressure  shock  and  a slightly  lower  pressure 
shock  are  greater  than  they  are  at  lower  pressures,  (in  the  lower  limiting  case, 
all  shocks  approach  sonic  velocity  and  never  overtake.)  For  most  tank  ruptures, 
airblast  pressures  are  comparatively  low  near  the  tank  so  that  a shock  front  that 

(C.l)  Swisdak,  M.”" Explosion  Effects  and  Properties;  Part  1 - Explosion  Effects 
in  Air",  Naval  Surface  Weapons  Center  NSWC/WOL/TR  75-116,  to  be  published. 
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starts  out  non-spherical  remains  non-spherical. 

This  leads  up  to  an  explanation  of  the  remarkable  differences  in  airblast 
overpressures  noted  in  Figure  C.5.  Both  the  15,000-psi  vessels  ruptured  along 
the  equatorial  weld,  separating  into  two  halves.  By  observing  the  direction  in 
which  the  two  haves  were  propelled,  it  is  possible  to  tell  which  side  of  the 
vessel  opened  first  since  the  two  halves  are  pushed  away  from  the  opening  in 
reaction  to  the  first  escaping  Jet.  The  airblast  pressures  will  be  higher  in 
the  vicinity  of  this  Jet. 

In  Figure  C.5,  Jet  directions  are  shown  in  relation  to  the  gage  lines  in  the 
upper  right  hand  inset.  Note  that  for  shot  1 (the  first  vessel  rupture)  the  Jet 
is  directed  along  the  vicinity  of  line  2.  Therefore,  blast  pressures  along  line  2 
are  higher  than  along  line  1.  (Line  3 was  not  in  use  on  this  test).  For  shot  2, 
the  Jet  is  directed  almost  directly  away  from  line  1 where  the  lowest  pressures 
were  recorded.  Similar  shot-2  pressures  are  read  along  line  2 and  the  gage  on 
line  3 since  the  Jet  was  about  equi-distance  between  the  two  lines.  Note  also 
that  the  shot-1  Jet  has  about  the  same  relation  to  line  1 as  does  the  shot-2  Jet 
to  line  2.  Thus  shot-1,  line-1  airblast  pressures  are  similar  to  shot-2  line-2 
airblast  pressures  in  Figure  C.5. 

Pressure  histories  from  the  shot  1 vessel  are  shown  in  Figure  C.6.  Positive 
shock  wave  durations  are  shown  in  Figure  C.7  and  positive  shockwave  impulses  are 
shown  in  Figure  C.8.  While  impulse  in  a strong  function  of  shockwave  overpressure, 
shockwave  durations  are  similar  for  all  pressure  levels  at  the  same  distance. 

Thus,  positive  shockwave  durations  in  Figure  C.7  are  not  as  strongly  related  to 
Jet  direction  as  is  the  overpressure-related  positive  shockwave  impulse  in 
Figure  C.8. 

C. 3 Airblast  from  30,000-psi  Vessel  Runture 

Measured  airblast  data  from  the  30,000-psi  vessel  ruptures  are  tabulated  in 
Table  C.3.  Airblast  overpressures  versus  ground  range  are  shown  in  Figure  C.9. 

The  data  displayed  in  Figure  C.9  are  again  typical  of  pressure  vessel-rupture 
data.  Pressure  amplitudes  correlate  with  the  estimated  Jet  orientations  shown  in 
the  Figure  C.9  inset.  The  shot-3  Jet  was  directed  along  line  2 where  highest 
pressures  were  recorded.  The  lowest  pressure  was  recorded  directly  away  from  the 
jet  by  the  single  gage  on  line  3. 

On  shot  7,  the  jetting  occurred  between  lines  1 and  3 giving  about  equal 
pressures.  Lowest  pressures  were  again  recorded  away  from  the  jet  on  line  1. 

Note  in  Figure  C.9  that  the  orientation  of  line  2 and  the  shot-7  jet  is 
similar  to  the  orientation  between  line  1 and  the  shot-3  jet.  Again  the  airblast 
overpressures  along  these  lines  are  similar. 

Pressure-time  histories  for  the  shot-3  rupture  are  shown  in  Figure  C.10. 

Positive  shockwave  durations  and  positive  shockwave  impulses  are  shown  in  Figures 
C.ll  and  C.12  respectively. 
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C.U  Airblast  from  50.000-psl  Vessel  Rupture 

Airblast  data  from  the  50,000-psi  vessel  ruptures  are  tabulated  in  Table  C.4. 
Airblast  overpressures  are  given  in  Figure  C.13.  Again, the  asymmetries  displayed 
by  the  data  in  Figure  C.13  are  typical  of  vessel  rupture  data. 

Both  vessels  seemed  to  jet  between  gage  lines  1 and  2.  However,  as  described 
in  Section  B.3»  the  shot-4  vessel  ruptured  with  the  fracture  plane  parallel  to  the 
ground;  the  shot-5  vessel's  fracture  plane  was  perpendicular  to  the  ground.  This 
may  account  for  the  relatively  small  difference  between  pressures  on  lines  1 and 
2 for  shot  4 and  the  correspondingly  larger  differences  for  shot  5»  The  reasoning 
here  is  that  the  shot-4  jet  would  have  been  wide  in  the  plane  parallel  to  the 
ground  but  narrower  in  the  vertical  plane;  the  reverse  would  be  true  of  shot  5. 
This  picture  may  be  further  illustrated  for  the  gages  at  a 1-ft  ground  range  and 
the  free  field  gage  at  1.75  feet.  For  shot  4 the  1-ft  ground  range  gage  shows 
higher  pressure  on  line  1 than  line  2.  However  the  line-l/line-2  difference  is 
not  as  great  as  for  shot  5*  Two  things  bring  this  about.  The  fact  that  the  shot- 
4 vessel  opened  in  a nearly  horizontal  plane  would  broaden  the  Jet  to  decrease 
differences  in  pressure  for  the  two  lines.  However,  the  rupture  plane  was  tilted 
to  favor  gage  line  1 where  the  highest  pressure  was  measured.  For  shot  5,  the 
rupture  plane  was  nearly  vertical  and  therefore  narrowed  the  Jet  in  this  plane. 
This  fact,  plus  the  slight  tilt  of  the  jet  toward  gage  line  2 accentuated  the 
differences  in  pressures  along  lines  1 and  2 for  shot  5. 

The  free  field  gages  were  up  near  the  tanks  center.  Thus  the  line-1  and 
line-2  gages  for  shot  4 were  both  almost  in  the  Jet,  and  line  1 for  shot  5 was 
also  near  the  Jet  and  registered  only  slightly  lower  pressure  than  the  above  two 
gages.  Its  companion  gage  on  line  2 for  shot  5 was  well  out  of  the  narrow  Jet 
and  registered  the  lowest  pressure. 

Pressure  histories  for  shot  5 are  shown  in  Figure  C.l4.  Positive  shockwave 
durations  and  positive  shockwave  impulses  are  shown  in  Figures  C.15  and  C.l6  re- 
spectively. The  large  scatter  in  positive  impulse  in  Figure  C.l6  reverses  the 
trend  of  the  30,000-and  15,000-psi  bursts  where  the  overpressure  scatter  was 
greater  than  the  scatter  in  impulse.  In  the  case  of  the  50,000-psi  rupture, 
the  impulse  scatter  is  greater  than  the  overpressure  scatter. 

C. 5 Airblast  from  the  15,000-psi  Confined  Vessel  Test 

As  shown  in  Figure  A. 6,  blast  gages  outside  the  confining  structure  walls 
were  laid  along  three  lines.  Airblast  pressures  measured  along  these  lines  are 
shown  in  Figure  C.17  and  are  given  in  Table  C.5.  Pressures  measured  along  the 
line  coming  out  of  the  box  (line  2)  are  higher  than  along  the  other  two  lines 
but  no  higher  than  was  predicted  for  a free  field  burst.  Pressures  behind  the 
box  are  slightly  lower  than  those  measured  at  one  side. 

An  attenuation  factor  relative  to  the  pressure  in  front  of  the  structure 
was  calculated  and  is  shown  as  a function  of  ground  range  and  overpressure  in 
Figure  C.l8. 
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Overpressure  ratios  for  I^/L^  range  from  a high  of  20  to  1 in  Figure  C.l8 
at  2.5  ft  from  ground  zero  to  2.2  to  1 at  the  60— foot  gage  positions.  This  shows 
greater  shielding  at  positions  near  the  structure  walls  while  the  shockwave  is 
less  affected  at  greater  distances. 

Such  behavior  is  expected  from  the  behavior  of  similarly  confined  conventional 
explosions.  Reference  C.3  reports  I^/L^  (front  to  side)  ratios  ranging  from  3:1 
at  20  psi  down  to  about  2:1  at  1 psi.  Lg/Lj  (front  to  back)  ratios  are  slightly 
lower,  3.5-'l  in  the  20-psi  region  down  to  2.2:1  in  the  l-psi  region.  The  front 
above  refers  to  the  structure 'e open  side.  Positive  shockwave  durations  are  given 
in  Figure  C.19.  The  uncertainties  in  measuring  positive  duration  along  line  3 
ruled  out  their  being  included  in  the  data.  (See  the  pressure  histories  in 
Figure  C.20.)  Measured  positive  shock  wave  impulses  are  given  in  Figure  C.21. 

Pressure  histories  measured  on  the  inside  walls  of  the  confining  structure 
are  shown  in  Figure  C.22.  (See  Figure  A. 6 for  the  location  of  gage  positions  1 
through  6.)  Based  on  the  analysis  in  Annex  B,  the  vessel  first  opened  on  the 
side  next  to  gages  5 and  6.  What  is  probably  an  airshock  reaches  gage  5 at  0.3 
milliseconds.  This  is  immediately  followed  by  a rise  which  may  be  the  argon  front. 
Almost  immediately,  the  gage  goes  off  scale,  probably  from  a parted  cable.  Gage 
6 shows  no  air  shock  but  does  show  a rise  probably  due  to  the  argon  front  before 
it  goes  off  scale;  again,  probably  due  to  the  cable  parting. 

Gages  3 and  h show  an  airshock  arrival  at  about  1 millisecond  after  burst. 

This  is  followed  by  a higher  pressure  at  about  1.5  milliseconds,  probably 
reflected  from  the  vessel  front  which  is  by  now  in  motion.  The  argon  front 
arrives,  driving  position  h off  scale  when  the  cable  parts  at  about  1.8  milli- 
seconds. The  signal  cable  for  gage  3 parts  at  about  3.2  milliseconds.  Note  that 

gage  It  is  6 inches  further  away  from  the  side  of  the  vessel  that  first  opened  than 

is  gage  3. 

Gages  1 and  2 see  a nominal  h.6  psi  shockwave  at  about  1 millisecond  after 
burst.  What  we  assume  to  be  the  argon  gas  front  arrives  at  about  3 milliseconds. 
Since  gages  1 and  2 are  on  the  side  opposite  the  side  where  first  rupture  occurred, 
the  low  air  shock  pressure  and  the  late  arrival  of  the  argon  may  be  expected. 

Based  on  the  above,  wall  load  pressures  are  from  two  sources:  the  airshock 

ahead  of  the  expanding  argon  gas  and  the  momentum  of  the  argon  itself.  In  addition, 

the  asymmetric  situation  within  the  confining  structure  results  in  multiple  re- 
flections at  the  structure  walls.  While  some  of  the  reflections  are  off  the  walls, 
others  may  be  off  the  argon  front.  Therefore  a qualitative  analysis  of  the  entire 
pressure-time  history  is  speculative.  The  situation  here  is  similar  to  that  of  a 
confined  conventional  explosion.  That  is,  an  airblast  wave  reaches  the  walls  first 
and  undergoes  multiple  reflections.  This  is  followed  by  the  static  pressure  due  to 
the  gaseous  products.  The  static  pressure  reaches  a peak  and  then  decays  at  a rate 
determined  by  the  venting  time  of  the  enclosure.  This  picture  fits  the  overall 
pressure  history  measured  in  our  case. 

(C. 3)  Keenan,  W.  and  Tancreto,  J.  "Blast  Environment  from  Fully  and  Partially 
Vented  Explosions  in  Cubicles",  Minutes  of  the  Sixteenth  Explosives  Safety  Seminar, 
Page  1527  Volume  II,  2b-2 6 Sep  1971*. 
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The  gross  pressure-time  functions  shown  in  Figure  C.22  are  probably  adequate  to 
determine  confining  structure  wall  response. 
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NSWC/WOL/TR  75-87 


8 


NSWC/WOL/TR  75-87 


GAGE  DISTANCE 
POSITION  (FT) 


L 

0 


1 

1 


1 

2 


i 

! 

i 


TIME  FROM  BURST  IN  MILLISECONDS 


t 


FIG.  C.4  PRESSURE  TIME  HISTORIES  FROM  478 -GRAM  : 

PENTOLITE  SPHERE;  HOB  = 1.6  FEET  \ 

■i 

1 


GAGE  DISTANCE 
POSITION  (FT) 

263 


TIME  FROM  BURST  IN  MILLISECONDS 


NSWC/WOL/TR  76-87 


GAGE  DISTANCE 

POSITION  {FT) 


5 MILLISECONDS  PER  DIVISION 


2 MILLISECONDS  PER  DIVISION 


FIG.  C.4  (CONTINUED) 


PREDICTED  OVERPRESSURE  FROM 
RUPTURED  VESSEL 


FREE  FIELD  MEASUREMENT 


GROUND  RANGE,  FEET 


FIG.  C.5  AIRBLAST  OVERPRESSURE  FROM  15,000- 
PSI  VESSEL  RUPTURE 
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FIG.  C.9  AIRBLAST  OVERPRESSURE  FROM  30,000  - PSI 
VESSEL  RUPTURE 
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FIG.  C.12  POSITIVE  SHOCKWAVE  IMPULSE  FROM  30,000  - PSI  VESSEL  RUPTURE 
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FIG.  C.14  PRESSURE-TIME  HISTORIES  FROM  50,000  PSI  SPHERE;  SHOT  5 
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FIG  C.16  POSITIVE  SHOCKWAVE  IMPULSE  FROM  50,000  - PSI  VESSEL  RUPTURE 
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FIGURE  C.19  POSITIVE  SHOCKWAVE  DURATION  FROM  15,000-PSI  CONFINED  VESSEL  RUPTURE 
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ANNEX  D FRAGMENT  MEASUREMENTS 


D.  1 Fragment  Instrumentation 

The  fragment  instrumentation  system  is  described  in  Annex  A.  The  fragment 
velocity  screens  and  their  placement  was  an  improved  system  based  on  the  one 
described  in  reference  3-  Unfortunately,  all  test  vessels  burst  into  two  pieces. 
As  a result,  they  seldom  hit  a breakwire  system;  and  when  they  did,  the  counters 
failed  to  stop.  Thus,  almost  no  fragment  velocities  were  measured  with  the 
screens. 


However,  in  almost  every  case,  one  portion  of  the  vessel  impacted  the  concrete 
firing  pad  with  considerable  force.  This  produced  an  n-wave  on  the  pressure-time 
(p-t)  histories  recorded  at  the  airblast  gage  positions  1-1,  1-2,  1-3,  and  2-3 
(Figure  A. 5)*  The  fragment  velocity  could  then  be  determined  from  the  arrival 
time  of  the  n-wave  and  the  distance  the  fragment  traveled.  Table  D.l  presents  the 
fragment  velocity  measurements.  Figure  D.l  shows  positions  of  recovered  fragments. 

D.2  Fragments  from  15,000-psi  Vessels 


Figure  D.2  shows  the  recovered  sections  of  shots  1,  2 and  6.  Figure  D.2a  is 
the  bottom  portion  of  the  shot-1  vessel  found  665  feet  from  GZ.  The  flatter  part 
is  where  the  fragment  impacted  the  concrete  pad.  The  recovered  section  weighed 
52.L  pounds.  The  averaged  measured  fragment  velocity  for  this  half  (before 

impact  with  the  concrete  pad)  was  318  feet/second.  Figure  D.2b  is  the  bottom  half 
of  the  shot-2  vessel  found  552  feet  from  GZ.  Again  the  flattened  portion  of  the 
fragment  is  that  part  that  struck  the  concrete  pad.  This  half  weighed  52.7  pounds. 
The  average  fragment  velocity  before  impact  with  the  pad  was  32U  feet/second. 

Figure  D.2c  and  D.2d  are  the  fragment  halves  from  the  confined  vessel,  shot  6. 
Figure  D.2c  is  the  bottom  portion.  Figure  D.2d  is  the  top  portion  of  the  vessel. 
The  top  portion  impacted  the  rear  wall  of  the  steel  confining  structure.  Both 
halves  weighed  51-1  pounds.  No  fragment  velocities  were  measured  for  this  test. 


For  the  15,000  psi  vessels  J m.,*-  were  C urieccvered  halves  - one  each  from 
shots  1 and  2.  Figure  D.l  shows  the  estimated  direction  these  halves  traveled. 

A search  for  the  unrecovered  fragments  out  to  1600  feet  for  shot  1 and  2500  feet 

for  shot  2 was  carried  out.  Since  the  search  area  for  the  shot-1  fragment  was 

!.  clear  plowed  ground,  the  lost  fragment  should  have  been  found.  That  it  was  not, 

; suggests  that  it  traveled  further  than  this  distance.  The  shot-2  fragment  path 

; was  over  an  area  covered  with  2-to  3- foot  high  grass.  Even  our  extended  search  for 

this  fragment  may  have  overlooked  its  position.  The  estimated  elevation  angles  of 
the  trajectory  was  5^*°  for  shot  1 and  6i*°  for  shot  2.  Calculation  using  a drag 
; coefficient  of  0.165  and  an  initial  velocity  of  325  ft/sec  show  a maximum  height 

and  distance  for  such  fragments  would  be  6l0  feet  and  1680  feet  respectively  for 
shot  1.  For  shot  2 the  maximum  height  and  distance  for  the  unrecovered  vessel 
s fragment  was  750  feet  and  1389  feet  respectively. 
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D. 3 Fragments  from  30.000-psi  Vessels 

Figure  D.3  shows  the  recovered  portions  of  shots  3 and  7.  Figure  D.3a  is  the 
bottom  half  of  the  shot-3  vessel  found  inside  the  arena.  This  section  weighed 
1 U 7 • 5 pounds.  The  mean  fragment  velocity  for  this  half  was  353  feet/second. 
Figures  D.3b  and  c are  the  ton  and  bottom  halves  of  the  shot-7  vessel  resnectively . 
The  average  fragment  velocity  was  2jb  feet/second  for  the  bottom  portion  of  the 
vessel.  (The  top  portion  did  not  impact  the  concrete  pad).  These  pieces  weighed 
9i.5  and  108  pounds  for  the  top  and  bottom  sections  respectively. 

The  unrecovered  fragment  from  shot  3 had  an  estimated  elevation  angle  of  70°. 
Figure  D.l  shows  the  azimuth  direction  of  this  half  of  the  vessel.  A search  of 
this  area  out  to  600  feet  was  carried  out.  Using  an  initial  velocity  of  353  feet/ 
second  for  this  fragment  the  calculated  maximum  height  and  distance  is  8l8  feet 
and  1136  feet  respectively.  Therefore,  the  fragment  probably  landed  in  a wooded 
area  where  we  were  unable  to  find  it. 

D.L  Fragments  from  50,000-psi  Vessels 

Figure  D.L  shows  the  recovered  fragments  from  Shots  L and  5.  Figure  D.La  is 
the  large  portion  of  the  shot-L  vessel  found  inside  the  arena.  This  section  of 
fragment  weighed  258  pounds  or  6 5%  of  total  tank  weight.  The  average  velocity 

from  this  fragment  was  216  feet/second  before  impact  with  the  concrete  firing  pad. 
Figure  D.Lb  and  c are  the  large  and  small  portions  of  the  shot-5  vessel.  Both 
pieces  hit  the  arena  wall.  The  large  piece.  Figure  D.3b,  weighed  271  pounds,  66 % 
of  the  total  weight.  The  small  portion  weighed  lL2  pounds.  The  larger  portion 
went  through  two  velocity  screens.  One  at  6.5  feet  indicated  a velocity  of  270 
feet/second.  The  other  screen  was  17.5  feet  away  and  was  tripped  at  6L.6  milli- 
seconds for  an  average  velocity  of  270  feet/second  also.  No  measurements  were 
made  on  the  smaller  portion. 

The  unrecovered  small  portion  of  the  shot-L  vessel  had  an  estimated  elevation 
angle  of  78°.  Figure  D.l  shows  the  direction  in  which  the  fragment  is  thought  to 
have  traveled.  A search  in  this  area  out  to  600  feet  was  carried  out.  Using  a 
drag  coefficient  of  0.165,  we  calculated  maximum  height  and  distance  of  886  feet 
and  722  feet  respectively. 
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FIG.  D.1  FIRING  SITE  AT  NSWC,  DAHLGREN,  VA. 
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ANNEX  E CALCULATIONS  OF  AIRBLAST  FROM  PRESSURIZED 

ARGON  SPHERES  ; 

(By  D.  Lehto) 

E.l  Equation  of  State  of  Argon 

Both  of  the  hydrocodes  used  here  (WUNDY  and  TUTTI)  require  an  equation  of 
state  either  in  the  general  form  P = P(p,E)  or  as  P(p,S)  isentropes  that  pass 
through  the  desired  initial  points.  The  Redlich-Kwong  equation  of  state 
(reference  E.l)  is  apparently  accurate  for  argon  (reference  E.2)  but  it  is  in  the 
P(p,T)  form  and  is  not  usable  in  our  hydrocodes  without  further  work.  If  a large 
number  of  calculations  involving  argon  were  to  be  done,  it  would  be  worth  while 
to  develop  a P(p,E)  equation  of  state.  Since  only  a few  initial  conditions  are 
of  interest  here,  we  take  the  simpler  course  and  use  individual  isentropes.  The 
expansion  isentrope  from  each  initial  condition  is  fitted  with  some  convenient 
function.  The  fits  used  here  are  given  in  Table  E.l. 

The  use  of  an  isentropic  argon  equation  of  state  (EOS)  is  not  quite  correct 
because  two  entropic  processes  occur  in  the  argon:  the  second  shock,  characteris- 

tic of  spherical  explosions,  that  moves  inward  from  the  expanding  gas/ambient  air 
interface,  and  any  reflected  shocks  from  nearby  walls. 

We  obtain  our  argon  EOS  data  from  the  tables  of  Din  (reference  E.3)  and  of 
Brahinsky  and  Neel  (reference  E.l*).  The  Din  tables  are  given  in  more  detail  in 
Little  (reference  E.5). 

Figure  E.l  shows  the  T-S  diagram  for  argon  with  lines  drawn  in  go  show 
isothermal  compression  to  15,000,  30,000  and  50,000  psi  and  the  isentropic 
expansion  paths  when  the  tank  is  ruptured.  When  the  argon  has  expanded  to  about 
1*0  atm  pressure,  the  saturation  line  is  reached  and  the  argon  will  become  opaque 
due  to  partial  condensation. 

Figure  E.2  shows  the  same  compression  and  expansion  paths,  this  time  on  a 
P-V  plot.  Note  that  the  isothermal  compression  curve  becomes  quite  steep  above 
about  15,000  psi  , indicating  that  little  energy  is  involved  in  further  compression 
to  50,000  psi  . This  means  that  a 50,000-psi,  17°C  tank  will  give  about  the  same 
airblast  as  a 15,000-psi  , 1J°C  tank.  The  expansion  isentropes  turn  sharply  at 
the  saturation  line. 


(E.l)  0.  Redlich  and  J.  N.  S.  Kwong,  "On  the  Thermodynamics  of  Solutions.  V.  An 
Equation  of  State.  Fugacities  of  Gaseous  Solutions,"  Physical  Review  Lh:  233-2hb 

(191*9). 

(E.2)  J.  P.  Kois,  personal  communication. 

(E.3)  F.  Din,  Thermodynamic  Functions  of  Gases,  Volume  2:  Air,  Acetylene,  Ethylene, 

Propane  and  Argon  ( Butterworth , London,  1956). 

(E.l*)  H.  S.  Brahinsky  and  C.  A.  Neel,  Tables  of  Equilibrium  Thermodynamic  Properties 
of  Argon.  Volume  III.  Constant  Entropy,  Arnold  Engineering  Development  Center 
AEDC-TR-69-19,  AD681*532,  Mar  1969. 

(E.5)  W.  J.  Little,  Tables  of  Thermodynamic  Properties  of  Argon  from  100  to  3000K, 
Arnold  Engineering  Development  Center  AEDC-TDR-6^-68,  AD!*3^*3d8,  Apr  196U. 
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One  of  the  main  purposes  of  these  calculations  was  to  provide  predictions  of 
blast  overpressures  for  the  experiments  to  be  done  with  15,000  and  50,000  psi 
tanks  at  ambient  temperature.  Consequently,  these  intial  conditions  are  stressed 

here. 

Figure  E.3  shows  the  P-V  isentropes  for  tanks  heated  to  1000°C  and  1750°C. 

These  curves  are  much  closer  to  ideal-gas  behavior  (straight  lines  of  slope  -5/3) 
than  are  those  for  the  cold  tanks  of  Figure  E.2.  The  point  for  50,000  psi  , 

1750°C  lies  in  a region  where  we  find  no  EOS  data.  Isentropes  in  this  region 
were  obtained  by  extrapolating  the  Brahinsky  and  Neel  (reference  E . U ) tables. 

This  extrapolation  is  not  very  accurate  because  of  the  strong  real-gas  effects  in 
this  region.  A Redlich-Kwong  EOS  adapted  to  P(p,E)  or  P(p,S)  would  be  useful  here. 

E.2  One-Dimensional  Flow  Calculations 

Figure  E.U  gives  the  calculated  peak  air  shock  overpressure  versus  distance 
from  the  tank  center  for  several  argon  pressures  and  temperatures.  The  solid 
parts  of  the  curves  are  "extrapolated"  pressures  from  the  WUNDY  (reference  E.6) 
calculations.  (An  artificial-viscosity  ca?culation  rounds  off  the  peak  of  a shock 
front.  The  true  peak  can  be  found  by  extrapolating  the  pressure  versus  radius 
profile  at  a fixed  time  to  the  radius  of  maximum  artificial  viscosity.  The  size 
of  the  extrapolation  is  usually  about  ten  percent  for  a WUNDY  calculation  and  may 
be  as  large  as  a factor  of  two  for  a TUTTI  calculation.) 

Within  about  one  tank  radius  from  the  tank  surface,  the  pressure  in  the 
expanding  argon  following  the  air  shock  can  have  a higher  pressure  than  the  air 
shock.  The  peak  argon  pressure  is  shown  by  dotted  lines  in  Figure  E.U. 

The  dashed  low-pressure  portions  of  the  curves  of  Figure  E. U are  obtained  by 
matching  the  TNT  free-air  curve  to  the  ends  of  the  solid  curves. 

The  points  at  a radius  of  0.62  feet  (the  initial  argon  sphere  radius)  are 
calculated  from  the  shock  tube  equation  as  shown  on  page  8l  in  reference  E.7. 

These  points  are  Joined  to  the  ends  of  the  solid  curves  with  smooth  dashed  curves, 
using  WUNDY  unextrapolated  pressures  (not  shown)  as  a guide. 

The  procedure  for  getting  the  shock-tube  initial  pressures  in  the  air  shock 
is  now  described  because  the  results  give  a useful  guide  in  interpolating  the 
: curves  of  Figure  E.U  to  pressures  and  temperatures  that  are  not  shown.  The 

procedure  is  as  follows: 

(1)  An  isentrope  is  prepared  for  the  desired  initial  argon  conditions. 

(2)  An  integration  is  carried  out  to  get  particle  velocity  versus  pressure 
in  the  isentropic  expansion. 

(e.£>)  D.  Eehto  and  M.  Lutzky,  One-Dimensional  Hydrodynamic  Code  for  Nuclear- 
Explosion  Calculations.  Naval  Ordnance  Laboratory  N0LTR  62-168,  AD615801,  Mar  1965. 
(The  WUNDY  Fortran  listings  in  this  report  are  obsolete.  This  report  is  being 
rewritten. ) 

(E.7)  W.  Liepmann  and  A.  Roshko,  Elements  of  Oas  Dynamics,  John  Wiley  and  Sous, 
Inc.,  London,  1958. 
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(3)  The  particle  velocity  versus  pressure  results  are  plotted  on  a graph  along 
with  the  shock  particle  velocity  versus  pressure  relation  for  the  desired  ambient 
gas  (here,  ideal  air),  as  shown  in  Figure  E.5.  The  intersection  of  these  curves 
gives  the  desired  air  shock  pressure  to  be  plotted  at  a radius  of  0.62  feet  on 
Figure  E.l*.  (if  both  the  argon  and  ambient  air  are  regarded  as  ideal  gases,  the 
air  shook  pressure  may  be  obtained  directly  by  solving  the  shock  tube  equation  by 
iteration. ) 

Figure  E.6  shows  this  initial  air  shock  pressure  varies  with  argon  pressure 
and  temperature.  At  high  temperatures,  the  argon  behaves  as  an  ideal  gas. 

Figure  E.7  shows  calculated  pressure  versus  time  calculations  at  distances  of 
3,  4,  5 and  6 feet  from  a one  cubic  foot  15,000  psia,  17°C  argon  tank.  The 
notation  "CS"  indicates  arrival  of  the  argon/air  contact  surface  at  the  gage 
location.  The  calculated  curves  are  somewhat  ragged  and  can  be  improved  by  using 
finer  zoning  in  the  WUNDY  calculation. 

E.  3 Two-Dimensional  Flow  Calculations 


The  results  described  so  far  are  for  a free-air  explosion  with  strictly 
radial  flow.  When  a reflecting  surface  is  present,  the  flow  is  no  longer  one- 
dimensional and  a two-dimensional  calculation  must  be  done.  A two-dimensional 
flow  Eulerian  hydrocode,  TUTTI,  was  written  for  this  purpose.  It  is  based  on  the 
FLIC  (Fluid-In-Cell)  method  of  Gentry,  et  al  (reference  E.8).  FLIC  as  described 
in  reference  E.8  is  for  only  one  material.  Air  and  argon  in  the  range  of  interest 
here  have  quite  different  properties  where  they  adjoin  at  the  contact  surface. 

The  argon  has  undergone  expansion  and  is  qui  -e  cold  (for  argon  initially  at  17°C 
in  the  tank),  while  the  air  has  undergone  compression  and  is  warm.  Continual 
computational  diffusion  of  these  two  materials  into  each  other  causes  error.  For 
this  reason,  TUTTI  was  written  as  a two-material  hydrocode  in  which  the  argon  and 
air  do  not  mix. 

We  now  consider  the  problem  of  a 15,000  psi  , one  cubic  foot  sphere  of  argon 
centered  5C  cm  above  a rigid  plane  surface.  The  argon  is  the  real  gas  of  Figures 

E.l,  E.2  and  E.3;  the  air  is  an  ideal  gas  with  a gamma  of  1.1*  (the  air  does  not 

get  compressed  enough  to  require  a real-gas  equation  of  state).  We  consider 
initial  argon  temperatures  of  17°C  and  1750°C. 

Figure  E.8  shows  TUTTI  results  for  the  shapes  of  the  contact  surface  and 
shocks  for  the  15*000- psi  , 17°C  tank.  The  reflected  shock  transmitted  into  the 
argon  is  not  shown;  it  lies  very  near  the  reflecting  surface  and  has  trouble 
moving  into  the  argon  because  of  the  strong  outward  argon  flow  and  because  of  the 
low  sound  speed  in  the  cold  mixed-phase  argon  (about  100  m/s).  The  reflected 
shock  tends  to  diffract  around  the  argon  rather  than  try  to  go  through  it.  This 

is  contrary  to  our  experience  with  H3  and  nuclear  explosions,  in  which  the 

reflected  shock  travels  more  rapidly  In  the  fireball  than  in  the  adjacent  air. 


(E.8)  R.  A.  Gentry,  R.  E.  Martin  and  E.  J.  Daly,  "An  Eulerian  Differencing  Method 
for  Unsteady  Compressible  Flow  Problems,"  J . Comput ■ Phys . 1 : 87-118  (i960). 
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Figures  E.9  and  E.10  show  the  TUTTI  pressure-time  calculations  along  the 
reflecting  surface.  The  large  rounded  peak  following  the  air  shock  for  the  17°C 
tank  in  Figure  E.9  is  due  to  the  pressure  in  the  expanding  argon.  This  peak  is 
represented  by  the  dotted  lines  in  Figure  E.^.  The  peak  air  shock  pressures  for 
the  17°C  tank  have  been  corrected  (dashed  curves)  by  applying  the  shock  reflection 
factors  of  reference  E.9  to  the  WUNDY  free-air  pressure-distance  curve  of  Figure 
E.L.  Figure  E.ll  shows  the  peak  air  shock  pressure  versus  distance  along  the 
surface.  The  rounding  off  of  the  TUTTI  air  shock  peak  is  considerable  but  lessens 
as  the  distance  from  the  tank  increases,  because  the  shock  grows  in  length  and 
thus  contains  more  of  the  calculational  zones.  Figures  E.12  and  E.13  show  the 
distribution  of  pressure  along  the  reflecting  surface  at  fixed  times. 

These  calculations  all  neglect  the  kinetic  energy  of  the  tank  fragments  and 
are  thus  expected  to  overestimate  the  shock  pressures  for  a spherical  tank 
rupture.  However,  the  actual  tanks  were  found  to  burst  open  like  clam  shells 
rather  than  break  up  into  small  fragments.  This  would  cause  the  pressures  to 
be  enhanced  in  the  direction  of  the  argon  Jet  from  the  splitting  tank  and  to  be 
reduced  in  the  opposited  direction.  The  Jet  is  expected  to  be  a broad  plume  in 
the  tank.  A rough  calculation  with  TUTTI  of  such  a plume  suggests  that  the 
shock  overpressures  in  the  direction  of  this  plume  will  be  about  a factor  of  two 
greater  than  in  a direction  90  degrees  from  the  plume  axis. 

These  TUTTI  calculations  can  be  made  arbitrarily  accurate  by  spending  an 
arbitrarily  large  amount  of  money  on  them,  but  the  cost  is  proportional  to  the 
cube  of  the  spatial  resolution.  The  TUTTI  runs  described  here  cost  about  $100 
each  and  occupied  130K  of  storage  on  a CDC  6500.  A 25  by  50  computational  grid 
was  used.  Because  the  reflecting  surface  was  a simple  plane,  it  was  possible  to 
sharpen  up  the  rounded  peaks  of  the  shocks  with  reflection  factors  applied  t ■>  the 
free-air  results.  For  more  complicated  shock  interactions  in  which  accurate  peak 
pressures  are  important,  more  accurate  TUTTI  calculations  are  needed.  If  it  is 
impulse  that  is  important  rather  than  peak  pressure,  the  rounding  of  the  shock 
does  little  harm. 
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IABL2  e.1  APJR0HMA1E  ANALYTIC  FITS  TO  THE  ARGON  ISENTROJES 

The  equations  are:  P = CQ(»r 

where  P a absolute  p re  a sure  (dyne/ca^) 
f ■ density  (g/cm’) ; 

log  P = Ci  ♦ C2(log  V)  * C3(log  V)2 

where  P = absolute  pressure  (psla) 

V ■ specific  volume  (cnP/mole) 
log  = log1Q. 
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FiG.  E.5  GRAPHICAL  SOLUTION  OF  SHOCK  - TUBE  EQUATION 
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FIG.  E.6 
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FIG.  E 7 FREE  AIR  PRESSURE  VS  TIME  AT  VARIOUS  DISTANCES  FROM  SPHERICAL  1 FTJ  15.000  PSI,  17°  C ARGON  TANK 
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FIG.  E 8 CONTACT  SURFACE  AND  AIR  SHOCK  CONFIGURATIONS  FOR  SPHERE  OF 
15,000 -PSI,  17°  C ARGON 
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FIG.  E.10  PRESSURE  VS  TIME  ALONG  THE  SURFACE  FROM  A 1 FT3  TANK  OF  ARGON  BURST  AT  15,000  PSI 
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FIG.  E-11  PEAK  AIR-SHOCK  OVERPRESSURE  VS  GROUND  RANGE  FOR  1 FT3  15,000  - PSI,  17°  C ARGON 
TANK  WITH  CENTER  50  CM  FROM  SURFACE 
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FIG.  E.12  PRESSURE  DISTRIBUTION  ALONG  SURFACE  FROM  SPHERICAL 
1 FT3  15,000  • PSI  ARGON  TANK  AT  17°  C 
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FIG.  E.13  PRESSURE  DISTRIBUTION  ALONG  SURFACE  FROM  SPHERICAL 
1 FT3  15,000  ■ PSIA  ARGON  TANK  AT  1750°  C 
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